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ABSTRACT 

Magnetic precursors of C-shocks accelerate, compress and heat molecular ions, modifying the kinemat- 
ics and the physical conditions of the ion fluid with respect to the neutral one. Electron densities are 
also expected to be significantly enhanced in shock precursors. In this Letter, we present observations 
of strongly polar ion and neutral molecules such as SiO, H 13 CO + , HN 13 C and H 13 CN, which reveal 
the electron density enhancements associated with the precursor of the young L1448-mm outflow. 
While in the ambient gas the excitation of the ions and neutrals is explained by collisional excitation 
by H2 with a single density of ~10 5 cm -3 , H 13 CO + shows an over excitation in the shock precursor 
component that requires H2 densities of a factor of >10 larger than those derived from the neutral 
species. This over excitation in H 13 CO + can be explained if we consider an additional excitation by 
collisions with electrons and an electron density enhancement in the precursor stage by a factor of 
^500, i.e. a fractional ionization of 5xl0 -5 . These results show that multiline observations can be 
used to study the evolution of the ion and electron fluids at the first stages of the C-shock interaction. 

Subject headings: stars: formation — ISM: individual (L1448) — ISM: jets and outflows — ISM: 
molecules 



1. INTRODUCTION 

In dense molecular clouds, the fractional ionization of 
the g as is very low (<1 0~ 7 ; Guelin, Langer & Wilson 
1982: iCaselli et al.lll998ft . Elect ron densities are there- 
fore insufficient to excite the first rotational levels of even 
strongly polar molecules like HCO + , HCN or HNC whose 
collisional cross sections (^10~ 6 -10~ 5 cm 3 s _1 ; Bhat- 



tacha r yya, Bhattacharyya fc N arayan 1981: ISaha et al.l 
IT98T1 iFaure Kr. TeuTv^Tnl Efiol are a factor of >100 
larger than th ose of weakly p olar species like CO 
(~1CT 8 cm 3 s" 1 : ISaha et al.lll981f) . All attempts to use 
molecular excitation to measure the effects of collisional 
excitation by electrons in t hese regions have been so far 
unsuccessful llLangerlll985j) . 

Modeling of the early stages of C-shocks predicts that 
the ion and elec tron densitie s are enhanced by the mag- 
netic precursor lpraindH9 80). For molecules like HCO+, 
HCN or HNC, the electron density enhancement by 
a factor of ~100 produced by the precursor (see e.g. 
iFlower et al.lll99ll) . would make electron collisions com- 
petitive with excitation by H2 collisions for the typical 
densities of dark clouds (~10 5 cm -3 ). Since the elec- 
tron collisional coefficients for high initial J transitions 
of HCO + are a factor of >10 larger than those of neu- 
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tral molecules like HC N at low temperatures (^10- 20 K; 
IChu k Dakarnol 119741 iBhattacharvva et aljll981l^ . dif- 
ferences in the molecular excitation between ion and neu- 
tral species are therefore expected in the precursor stage 
6 . Toward the young L1448-mm molecular outflow, the 
detection of very narrow SiO emission and the enhance- 
ment of the ion abundance have been interpreted as 
signa tures of the shock precursor l)Jimenez-Serra et alJ 
2004). Multiline observations of strongly polar species 
like HCO+, HNC and HCN toward this outflow, are ex- 
pected to show differences in excitation between the am- 
bient gas where electron excitation is negligible, and the 
shock precursor where electron collisions should be im- 
portant. 

In this Letter, we present observations of several rota- 
tional transitions of SiO, H 13 CO+, HN 13 C and H 13 CN 
observed toward the ambient and shock precursor com- 
ponents of the L1448-mm outflow. The excitation differ- 
ences observed between H 13 CO + , and SiO, HN 13 C and 
H 13 CN in the ambient and precursor components, can be 
explained by the electron density enhancement expected 
at the first stages of the C-shock evolution. 

2. OBSERVATIONS & RESULTS 

We observed several transitions (from J —1 to 5) of 
SiO, H 13 CO+, HN 13 C and H 13 CN toward three different 
positions in the L1448-mm outflow. Except the J = 4^3 
lines of H 13 CO+, HN 13 C and H 13 CN, all the molecular 
transitions were observed with the IRAM 30 m telescope 
at Pico Veleta (Spain). We used the wobbler-switched 
and frequency-switched modes with position and fre- 
quency throws of 240" and 7.2 MHz. The beam size was 
~28", 19" and 11" at -90, 130 and 260 GHz. The SIS 

6 These excitation differences between ion and neutral molecules 
reflect a fundamental difference between the corresponding 
electron— impact cross sections at low energy: for neutrals, cross 
sections g o to zero at threshold w hile for ions, they are large and 
finite (e.g. [Chu & Dalgarno 1971). 
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Fig. 1— Low- and high-J lines of SiO, H 13 CO+, HN 13 C 
and H 13 CN observed toward L1448-mm. Offsets in arcseconds arc 
shown in the upper part of the columns (oj(j 2 ooo) = 03 h 25 m 38 s 0, 
5(j2ooo) = 30°44'05"). The vertical dotted lines show the ambi- 
ent gas at 4.7kms — 1 and the precursor component at 5.2kms~ 1 
toward L1448-mm (0,0), and at 5.4 and 5.6kms _1 toward L1448- 
mm (0,-10) and (0,-20), respectively. 



receivers were tuned to single-side band with rejections 
of >10dB. We used the VESPA spectrometers with a 
spectral resolution of 40 kHz, i.e. a velocity resolution of 
-0.14, 0.09 and 0.05 km s" 1 at -90, 130 and 260 GHz. 
Typical system temperatures ranged from 100 to 500 K. 

The J = 4^3 lines (-347 GHz) of H 13 CO+, HN 13 C 
and H 13 CN were observed with the JCMT telescope at 
Mauna Kea (Hawaii) in the frequency-switched mode 
with a frequency throw of 16 MHz. The beam size 
was —14", which matches the IRAM 30 m beam for the 
J =3^2 lines of H 13 CO+ and HN 13 C. We used the B3 re- 
ceiver in dual-mixer and single-side band mode with an 
image rejection of 12-14 dB. The DAS spectrometer pro- 
vided a spectral resolution of 156 KHz (~0.14kms _1 ). 
The system temperatures were of 450-560 K. All the in- 
tensities were calibrated in units of antenna temperature 
and converted to main beam temperatures using efficien- 
cies of 0.82, 0.74 and 0.52 at -90, 130 and 260 GHz for 
the IRAM 30 m data, and 0.63 for the JCMT data. 



Fig. 1 shows the line profiles of all transitions mea- 
sured toward L1448-mm (0,0), (0,-10) and (0,-20), and 
Tab. 1 gives the observed parameters for the different 
velocity components. As expected for a shock tracer 
in the precursor (Jimenez-Serra et al. 2004), the nar- 
row (— 0.6kms _1 ) SiO J — 2— >1 and 3^2 lines have 
single-peaked profiles whose peak emission is slightly 
rcdshiftcd (— 5.2kms ; Fig. 1) with respect to the am- 
bient 4.7kms- 1 cloud. The J = 1^0 lines of H 13 CO+, 
HN 13 C and H 13 CN show double-gaussian profiles (the 
ambient and shock precursor components) and have 
linewidths of — 0.6-0. 7 km s -1 for each of the velocity 
components. The H 13 CO + emission peaks at the shock 
precursor component toward the positions where narrow 
SiO is detected (Jimenez-Serra et al. 2004). However, 
the HN 13 C peak emission is centered at the ambient 
cloud in L1448-mm (0,0). Toward L1448-mm (0,-10) 
and (0,-20), the HN 13 C emission is progressively red- 
shifted to 5.4 and 5.6kms _1 , respectively (see vertical 
doted lines in Fig. 1). H 13 CO+ and H 13 CN also peak at 
5.4 and 5.6 km s^ 1 toward these positions, as if we were 
observing the evolutionary effects of the propagation of 
C-shocks through the unperturbed gas. The detection of 
broad SiO emission (terminal velocity of —25 kms -1 ) to- 
ward L1448-mm (0,-10) and (0,-20), supports this idea. 

The J = 3->2 and 4^3 lines of H 13 CO+, HN 13 C and 
H 13 CN tend to be single-peaked and have linewidths 
of -O^kms- 1 (Tab. 1). The HN 13 C emission arising 
from the ambient gas toward L1448-mm (0,0) is weak. 
However, the high-J H 13 CO + emission measured to- 
ward L1448-mm (0,0) and (0,-10), is very bright and 
mainly arises from the precursor component. H 13 CO + 
also shows faint emission centered at the ambient cloud 
toward L1448-mm (0,0). The SiO J = 5->4 and the 
H 13 CN J = 4— >3 lines have not been detected in any of 
the velocity components toward L1448-mm. 

3. EXCITATION DIFFERENCES BETWEEN ION AND 
NEUTRAL MOLECULAR SPECIES 

A first look to the high-J lines of H 13 CO+, HN 13 C 
and H 13 CN, clearly shows that the emission of H 13 CO + 
is anomalously bright in the shock precursor component 
compared to that of the neutrals. Since the emission of 
all these species is expected to be optically thin, the line 
intensity ratio between different transitions is directly re- 
lated to the excitation temperature of the ion and neutral 
fluids in the ambient and shock precursor gas. The in- 
tegrated line ratios between the J = 3— >2 and 1— >0 lines 
for H 13 CO+ and HN 13 C, and between the J = 4^3 and 
l->0 lines for H 13 CO+ and H 13 CN in the ambient and 
precursor components, are shown in Tab. 2. In the ambi- 
ent cloud, the line ratios of H 13 CO + are very similar to 
those of HN 13 C and H 13 CN. However, the line intensity 
ratios of H 13 CO + in the precursor component are up to 
a factor of 9 larger than those of the neutrals, indicat- 
ing a higher excitation for the ions in this component. 
H 13 CO + is "selectively" excited in the precursor gas. 

We can estimate the H2 densities required to explain 
the line intensity ratios of Tab. 2, by using a model for 
the excitation of the observed molecules. Let us con- 
sider the LVG approximation and the only excitation by 
H? collisions . We have use d the H? collis ional rates of 
iTurner et al] lH99l for S iO.lFloweJ 111991 for H 13 CO+, 
and IGreen fc Thaddeusl (I1974D for HN 13 C and H 13 CN. 
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TABLE 1 







(0,0) 






(0,-10) 






(0,-20) 




Line 


V LSR 


Av 


T A * 


V LSR 


Ad 




V LSR 


Av 


T A * 




(kms -1 ) 


(kms -1 ) 


(K) 


(kms -1 ) 


(km s — ) 


(K) 


(kms" 1 ) 


(km s~ ) 


(K) 


Si0(2->1) 


—4.7 




<0.012 


—4.7 




<0.048 


—4.7 




<0.060 


5.170(8) 


0.62(2) 


0.107(5) 


5.18(3) 


0.60(8) 


0.11(2) 


5.11(2) 


0.41(8) 


0.10(2) 


SiO(3->2) 


-4.7 




<0.084 


-4.7 




<0.132 


-4.7 




<0.183 


5.08(7) 


0.6(2) 


0.14(1) 


5.02(6) 


0.5(1) 


0.17(3) 


5.20(7) 


0.9(1) 


0.18(4) 


SiO(5^4) 


-4.7 




<0.171 


-4.7 




<0.130 


-4.7 




<0.084 


-5.2 




<0.171 


-5.2 




<0.130 


-5.2 




<0.084 


H 13 CO+(1-+0) 


4.584(7) 


0.59(1) 


0.629(8) 


4.685(6) 


0.70(1) 


0.62(2) 


4.84(2) 


0.84(5) 


0.46(3) 


H 13 CO+(3^2) 


5.300(5) 


0.81(1) 


1.160(8) 


5.424(0) 


0.796(8) 


1.26(2) 


5.545(9) 


0.84(2) 


1.00(3) 


4.7(0) 


0.6(0) 


0.28(8) 


-4.7 




<0.168 


-4.7 




<0.225 


H 13 CO+(4^3) 


5.356(6) 


0.69(2) 


2.25(8) 


5.48(2) 


0.95(5) 


0.66(8) 


-5.2 




<0.225 


4.7(0) 


0.55(0) 


0.08(5) 


-4.7 




<0.162 


-4.7 




<0.138 


HN 13 C(1^0) 


5.338(7) 


0.62(2) 


1.43(5) 


5.47(2) 


0.63(4) 


0.52(4) 


-5.2 




<0.138 


4.785(3) 


0.806(7) 


0.72(1) 


4.81(2) 


0.65(2) 


0.56(2) 


4.85(2) 


0.50(4) 


0.34(3) 


HN 13 C(3^2) 


5.468(6) 


0.78(2) 


0.36(1) 


5.53(2) 


0.88(3) 


0.71(2) 


5.59(1) 


0.86(4) 


0.68(3) 


4.75(5) 


0.6(1) 


0.35(4) 


-4.7 




<0.261 


-4.7 




<0.264 




-5.2 




<0.237 


-5.2 




<0.261 


-5.2 




<0.264 


HN 13 C(4^3) 


4.38(8) 


0.8(0) 


0.16(3) 


-4.7 




<0.117 


-4.7 




<0.108 




-5.2 




<0.144 


-5.2 




<0.117 


-5.2 




<0.108 


H 13 CN(1^0) 


4.86(5) 


1.06(7) 


0.18(2) 


4.89(2) 


0.80(6) 


0.27(1) 


4.9(1) 


0.6(2) 


0.16(1) 


H 13 CN(4->3) 


5.32(2) 


1.25(4) 


0.21(2) 


5.52(2) 


0.70(6) 


0.30(1) 


5.5(1) 


0.7(2) 


0.28(1) 


-4.7 




<0.141 


-4.7 




<0.138 


-4.7 




<0.141 


-5.2 




<0.141 


-5.2 




<0.138 


-5.2 




<0.141 



TABLE 2 
(0,0) (o-io) 





Amb. 


Pre. 


Amb. 


Pre. 


H 13 CO+ 


0.5 


1.6 


<0.1 


0.6 


HN 13 C 


0.4 


<0.2 


<0.2 


<0.1 


H"CO + 


0.1 


0.9 


<0.1 


0.3 


H 13 CN 


<0.2 


<0.1 


<0.2 


<0.2 



Considering that the emissions of all molecules have sim- 
ilar spati al distributions, a nd a kinetic temperature of 
21 K (see iCuriel et al.l lT999) . the estimated H 2 densities 
and molecular column densities are shown in Tab. 3. For 
the ambient gas, the derived H2 densities for all molecules 
are of few 10 5 cm -3 , consistent with excitation by H2 col- 
lisions. The H2 densities derived from SiO, HN 13 C and 
H 13 CN in the precursor are similar to those of the am- 
bient gas. However, as expected from the large H 13 CO + 
line ratios, the H2 densities required for this ion in the 
precursor gas are a factor of >10 larger than those for the 
neutral species (Tab. 3). This clearly illustrates that ex- 
citation only by H2 collisions with a single density cannot 
explain the excitation of H 13 CO + in the precursor. 

4. COLLISIONAL EXCITATION BY ELECTRONS. 

Since the ions have been selectively excited in the pre- 
cursor component by an extra mechanism beside the H2 
impact, we explore the possibility that this selective exci- 
tation is produced by collisions with electrons. The effi- 
ciency of excitation of molecular ions by electrons can be 
significantly larger than that of neutral molecules at the 
low temperatures of dark clouds. To illustrate this, we 
compare the electron collisional rates of HCO + and HCN 
for the J = 0— >1 and 1^2 transitions at low and high 
temperatures. While the HCO + /HCN collisional coeffi- 
cient ratio is only -4.6 at 100 K, this ratio is increased to 
~6 at a temperature of 10 K ijBhattacharvva et al.lll98ll 



iSaha et al.lll981|) . This difference between the HCO+ 
and HCN rates is expected to furth er increase for higher 
initial J and large AJ transitions ifBhattach arvva et alJ 
[MU). New calculations of the HCO+ and H 13 CO+ col- 
lisional rates for all trans itions between J = 1 and J = 5 
l|Faure fc Tennvsonl I2006T) . show that the HCO+ rates 
with AJ>2 exceed those of HCN by more than one order 
of magnitude at 10 K. This naturally introduce a differen- 
tial excitation between the ions and neutrals as observed 
in the precursor component. 

We can constrain the electron density required to re- 
produce the H 13 CO + line intensities observed in the pre- 
cursor component by using the LVG model including col- 
lisions with both H2 and electrons. For an H2 density 
of -^3xl0 5 cm~ 3 (similar to that derived from HN 13 C 
and H 13 CN; Tab 3) and a temperature of 21 K, the es- 
timated electron densities in the precursor component 
toward L1448-mm (0,0) and (0,-10) in the optically thin 
case [N(H 13 CO+)-10 n cm- 2 ], are ~240 and 600c m - 3 
respectively, which correspond to fractional ionizations of 
8x 10 -4 and 2 x 10 -3 . However, in the optically thick case 
[N(H 13 CO+)~7-9xl0 12 cm- 2 ], and for higher tempera- 
tures (-^35-45 K), the fractional ionization in the precur- 
sor decreases to -^5 x 10~ 5 . Although even higher temper- 
atures (~100K) could reproduce the H 13 CO + intensities 
in the precursor component for a fractional ionization of 
<10~ 7 , the derived HN 13 C and H 13 CN line intensities 
would clearly exceed (by up to a factor of 5) the upper 
limits of Tab. 1. The derived ionization fraction implies 
an electron density enhancement by a factor of ~500 with 
respect to that of the quiescent gas (<10~ 7 ). 

Cons idering an extrapo lation of the electron collisional 
rates of lSaha et alJ l)1981|) for HCN, we can now estimate 
the expected line intensities of HN 13 C and H 13 CN in the 
precursor for the fractional ionization (~5x 10~ 5 ) derived 
from H 13 CO + . The expected intensities are similar to 
those reported in Tab. 1, except for the HN 13 C J = 3^2 
line whose predicted intensity exceeds the upper limits 
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in Tab. 1 by a factor of ~2. Given the uncertainties in 
the rates, and the relative spatial distribution of the ion 
and neutral gas, the data are consistent with the idea of 
an electron density enhancement in the precursor. High 
angular resolution observations are required to establish 
the spatial distribution of the ion and neutral species in 
this component. 

5. ON THE ORIGIN OF THE ELECTRON DENSITY 
ENHANCEMENT 

Toward the quiescent gas of L1448-mm, the ion and 
neutral fluids show similar excitation conditions. In fact, 
the H2 densities obtained from the ions and neutrals 
are all consistent with few 10 5 cm -3 for this component. 
Since the fractional ioni zation is expected to be of <10~ 7 
in th e ambient cloud l)Guelin et alJ [1982: iCaselli et al.l 
1199^1 . the high- J H 13 CO+ and HN 13 C excitation in this 
component is completely dominated by H2 collisions. 

In contrast with the quiescent gas, the precursor com- 
ponent shows an over excitation in H 13 CO + . The line 
ratios and H2 densities estimated for this ion in the pre- 
cursor gas, are a factor of 10 larger than those for SiO, 
HN 13 C and H 13 CN. In section 4, we have shown that 
an electron density enhancement by a factor of ^500 
in this component could explain the over excitation in 
H 13 CO+. Modeling of C-shocks shows that the UV flu- 
orescence radiation generated by the collisional excita- 
tion of H2 in the magnetic precursor rapidly enhances 
the ion and electron densities in thi s region (t<100yrs 
from the inception of the C-shoc k; iFlower et alJ Il996t 
iFlower fc Pineau" des Foretsl 12003(1 . One may think 
that the probability to detect this enhancement toward 
molecular outflows for such short time-scales is negli- 
gible. However, L1448-mm is a very young outflow 
(tdi/n~1000yrs), and the dynamical time-scales derived 
from the proper mot ions of the SiO bullets (~90yrs; 
iGirart fc Acordll2001[) are consistent with the possibility 
of detecting the shock precursor as predicted by C-shock 
models. In fact, the more redshifted velocities toward 
L1448-mm (0,-10) and (0,-20) are consistent with the 
observation of different evolutionary stages of C-shocks 
in the different positions in the outflow as predicted by 
models. 

We can use our results to constrain the electron den- 
sity enhancement produced by the precursor. For op- 
tically thin emission and low temperatures, the derived 
fractional ionization is ~10~ 3 . Flower et al. (1996) and 
Flower & Pineau des Forets (2003) predicted that the 
ionization fraction is increased to ~10~ 5 in the precursor 
stage. Our estimate of the fractional ionization clearly 
exceeds these results by a factor of 100, and even ex- 
ceeds the cosmic abundance of atomic carbon (the main 
repositor y of positive charge in dark clouds; x(C)^2- 
3xl0~ 4 ; ICardelli et all 11996ft by a factor of -5. How- 
ever, if we increase the kinetic temperature (i ons are 
expected to be rapidly heated by the precursor; iDrainei 
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1980) and consider optically thick emission, the ioniza- 
tion fraction can be decreased to ~5xl0~ 5 . This result 
is consistent with the model predictio ns l|Flower et alJ 
ll996tlFlower fc Pineau des Foretsll2003|) . 

We cannot rule out the possibility that the ion 
and electron enhancement in the precursor component 
is produced by the radiative precursor of J-shocks 
ijShull fc McKeeill979|) . Chemical models that include 

TABLE 3 



Molecule H-2 density (cm ^) Column density (cm 







(0,0) 






Ambient 


Precursor 


Ambient 


Precursor 


H 13 CO+ 


1.3xl0 5 


2.0xl0 6 


6.7X10 11 


2.2X10 12 


SiO 




2.7xl0 5 




1.4x10" 


HN 13 C 


1.9xl0 5 


<3.4xlO s 


1.4xl0 12 


5.1x10" 


H 13 CN 


<1.8xl0 6 


<1.6xl0 6 


2.1xl0 n 


2.4x10" 






(0,-10) 




H 13 CO+ 


<1.3xl0 5 


6.0x10 s 


4.2x10" 


1.8xl0 12 


SiO 




4.0x10 s 




1.3x10" 


HN 13 C 


<1.9xl0 5 


<1. 1x10 s 


l.lxlO 12 


1.9xl0 12 


H 13 CN 


<1.3xl0 6 


<1. 2x10 s 


3.0x10" 


3.4x10" 


H 13 CO+ 




(0,-20) 


l.OxlO 12 


<1.6xl0 5 


<5.7xl0 4 


3.1x10" 


SiO 




4.5x10 s 




1.3x10" 


HN 13 C 


<4.0xl0 5 


< 1.3 x10 s 


4.5x10" 


1.6xl0 12 


H 13 CN 


<2.0xl0 6 


<1. 2x10 s 


1.9x10" 


3.2x10" 



illumination by UV photons predict the enhancement of 
molecules such as HCO+, HCO and H CN for t<300yrs 
(IViti fc WilliamsllT999tlViti et al]l200l How ever. HCO 
(a typical PDR tracer; iSchenewerk et al.lll988fl is not de- 
tecte d in the precursor component l|,Iimenez-Serra et alJ 
2004), which suggests that the ion and electron enhance- 
ments are likely due to the magnetic precursor of C- 
shocks. 

In summary, the differences in the kinematics and ex- 
citation between the ion and neutral components in the 
L1448-mm molecular outflow are clear indicators of the 
early interaction of C-shocks with the ambient gas. The 
over excitation in H 13 CO + has allowed to measure, for 
the first time, the electron density enhancement in the 
precursor of a C-shock. The estimated fractional ioniza- 
tion in the precursor component is of ^5xl0~ 5 , which 
implies an enhancement of the electron densities by a 
factor of >500 with respect to the ambient gas. 
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